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Abstract 

We confront the lattice results on QCD phase diagram for two and three 
flavors with the hadron resonance gas model. Taking into account the trunca- 
tions in the Taylor-expansion of energy density e done on the lattice at finite 
chemical potential fj,, we find that the hadron resonance gas model under the 
condition of constant e describes very well the lattice phase diagram. We also 
calculate the chemical freeze-out curve according to the entropy density s. The 
s-values are taken from lattice QCD simulations with two and three flavors. We 
find that this condition is excellent in reproducing the experimentally estimated 
parameters of the chemical freeze-out. 



1 Introduction 

At high temperatures T, it is conjectured that the hadrons are dissolved into quark- 
gluon plasma (QGP). Reducing the temperature, QGP will hadronize. At T c h, the 
system goes into chemical equilibration, freezing-out. At T < T c h, thermal equilibra- 
tion will take place and the matter consists, once again, of non-interacting hadrons. 

In this work, we confront the lattice results on QCD phase diagram with the 
hadron resonance gas model (HRGM) [H |2] . We study the influence of strange quarks 
on the location of the phase diagram. To compare with two flavor lattice results, 
we include in HRGM the non-strange resonances up to mass 2 GeV. For the three 
flavor results, we include all observed resonances with masses up to 2 GeV. In the 
Boltzmann limit, the particle number, energy and entropy density for one particle and 
its anti-particle, respectively, read 
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These thermodynamic quantities will be summed over all resonances taken into ac- 
count. In mapping out the phase diagram, we take advantage of our previous work [1 
on simulating T c for different quark masses and reproducing lattice thermodynamics 
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Tab. 1: QCD phase transitions from lattice QCD simulations with different flavors at \i = 0. 

at zero and finite \x. We assume that the critical energy density e c (T c ,/i = 0) remains 
constant at all /x-values. We suppose that this assumption is not affected by the 
existence of different transitions (Tab. [TJ. As done on the lattice, we can stop the 
Taylor-expansion of trigonometric function in Eq. (J2J) up to certain terms. Taking 
into account the quantum statistics, the results are drawn in Fig. ^ and Fig. El 

The lattice results at /i = are summarized in Tab. ^ We find that the strange 
quarks, due to the heavy mass, have small effect on the critical values. At \i ^ 0, 
lattice simulations still suffer from the sign-problem. The fermion determinant gets 
complex and therefore MC techniques are no longer applicable. However, considerable 
progress has been made to overcome this problem EJ |H| • 

For the chemical freeze-out curve, we compare its experimentally estimated pa- 
rameters, T and fib, with HRGM results. We propose that the entropy density s is 
the thermodynamic condition which drives the chemical freeze-out. Its value is taken 
from lattice calculations for different flavors at /i = 0. The resulting curve is seen 
in Fig. El Apparently, this condition gives excellent agreement with the experimen- 
tal data and reproduces very well the two characterizing endpoints of the freeze-out 
curve, (T = 0,/x ^ 0) and (T ^ 0,/i = 0). 



2 The QCD phase diagram 

In lattice calculations, the dependence of T c on /i can be determined from the pseudo- 
critical coupling /3 c (/i). T c is the value at which the susceptibility becomes maximum. 
From the first non-trivial Taylor coefficients of T c (fi) 

* ;TM = -"^m* . (4, 



dfi 2 cy ^' T c (ji = Q) <9/i 2 V dP 

N T and a are the temporal lattice dimension and the lattice spacing, respectively. The 
lattice beta function (3{a) is needed to get the physical units. The results for two jH] 
and three flavors jH], respectively, are 

1 0.070(35) (-—^ T Y. (5) 



T c {n q = 0) ■ v ' \T c (fi q = 0) 

= 1 - 0.114(46) ( — -] , (6) 

where \i q = ^1 '3 is the quark chemical potential. ^ is the baryo-chemical potential. 

In Fig. m we plot Eq. (0) as vertical lines. To compare with two flavor lattice 
results, we entirely exclude the strange hadron resonances. The curves give the results 
obtained from different truncations in the Taylor-expansion of e c , Eq. (J2J). The three 
flavor results are shown in Fig. El Here, we explicitly set the strange quark chemical 
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Fig. 1: T c vs. lib for light quarks. We com- 
pare lattice with HRGM results. The truncation 
up to the second-order is obviously able to de- 
scribe both lattice simulations |SJ|nj. The solid 
line gives our predictions for the phase diagram, 
when the e c -expansion is not truncated. 
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Fig. 2: The same as in Fig.^ Here, we include 
the strange quarks. Up to second-order trunca- 
tions of e c , Eq. @ gives a good description for 
the lattice results The region of most 

reliable lattice results is < ^ < 0.5 GeV. The 
solid circle is the critical endpoint |1U| . 



potential /i s = 0, as the case in the lattice calculations. An extensive discussion about 
the dependence of fi s on T and ^ is given in Ref. [2]. 

We also show the latest lattice calculations of the critical endpoint. According to 
Ref. PU, it is located at /if = 360 ± 40 MeV and = 162 ± 2 MeV. In Ref. P], it 
has been found that = 420 MeV. The T c ep -coordinate has not yet been calculated. 
We assume that the critical endpoint, at which the phase transition is second-order, 
does not affect our results. 



3 The chemical freeze-out 

As discussed in Ref. [12], the entropy density is the thermodynamic condition which 
guarantees chemical equilibration between reactants and products. To prove this, we 
recall the theoretical chemistry. Without energy input the chemical reactions always 
proceed toward equilibrium. The equilibrium constant C is related to the energy 
difference between reactants and products 5e via Boltzmann factor, C ~ exp(<5e/T). 
From second law of thermodynamics, 5e = 5F + T5s + fi5n. 5F is the free energy 
difference between reactants and products. It represents the total work in the reacting 
system. At equilibrium, 5F = and, therefore, 

C w exp -s-n (^J . (7) 

The entropy gives the amount of energy which can't be used to produce further work. 
Then, the equilibrium entropy reads 



s(T,p) « lnQ) — n(T, y) (£) 



In literature, there are different models E31 E] for the freeze-out conditions. 
Cleymans and Redlich ^H] assumed that C depends on some average energy per hadron 
particle. This apparently ignores the change in the baryon number with increasing 
/Lt. Braun-Munzinger and Stachel [H] rectified this and assumed that the freeze-out 
is given by a constant baryon number density (compare this with Eq. (|7|)). The two 
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Fig. 3: The freeze-out curve according to the entropy density, Eq. J2J. The two endpoints (large 
solid circles) and the experimentally estimated results (squares) are very well reproduced. There is 
no difference between strange and non-strange hadron resonances. The QCD phase diagrams for two 
and three quark flavors n / are also drawn. 

models in Ref . [T^l IT5] give identical results. As we show in Ref . jE| , there are some 
discrepancies in reproducing the experimental data by all these models. 

We calculate the freeze-out curve according to Eq. 0. In this context, the entropy 
can be seen as a measure for the degree of sharing and spreading the energy inside 
the system. The way of distributing the energy is not just an average value. But the 
method that controls the chemical equilibration, i.e., produces no additional work. 
This is the equilibrium entropy, Eq. (JEl). 

As T — > and \i ^ 0, it is assumed that the system is just a degenerate Fermi gas 
of nucleons. Then from Eq. (j3*jl. s — 0. At /x = and T ^ 0, the system becomes 
degenerate Bose gas of pions and rho-mesons. Therefore, s gets a finite-value. In the 
intermediate region, along the freeze-out curve, s ranges between these two extreme 
limits. 

At small /i, it is supposed that the freeze-out and phase transition are coincident. 
We can therefore use the lattice calculations [21 CI ; s/T 3 = 5 for = 2 and 
s/T 3 = 7 for rif = 3. The normalization with respect to T 3 should not be connected 
with the massless ideal gas. Either the hadrons in HRGM or quarks on lattice are 
massive. The constant ratio s/T 3 obviously fulfills Eq. (JHJ) and simultaneously obeys 
the third law of thermodynamics. The quantum entropy ^EJ is not included here. 
On the other hand, assignment s/T 3 to a constant value and numerically solving 
Eq. with respect to T, leads to the excellent results seen in Fig. El The two 
characterizing endpoints and experimental data are very well reproduced. There is 
almost no difference between including and excluding strange resonances, as long as 
we use the corresponding s-value. 

4 Conclusion 

We compared between the lattice results on the QCD phase diagram for different 
flavors and the HRGM results. Taking into account the truncations in the Taylor- 
expansion of the energy density e done on the lattice, we found the condition of 
constant energy density results in an excellent agreement with all available lattice 
results. We conclude that the influence of the strange quarks on the location of QCD 



phase diagram is relatively large at small baryo-chemical potential At large //&, 
there is almost no influence. This conclusion is valid unter the assumption, that the 
strange quark chemical potential fi s is vanishing. 

For mapping out the chemical freeze-out curve, we used the entropy density. Tak- 
ing its value from lattice simulations with two and three quark flavors at \i — and 
assuming it remains constant on the entire /z-axis, we obtain the meaningful results 
shown in Fig. |3J The experimentally estimated data is seen to be very well described 
under this condition. The two characterizing endpoints of the chemical freeze-out 
curve are also reproduced. 
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